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Abstract 

Fully  Automatic,  hp-Adaptive  Simulations  for  Maxwell’s  Equations 

Leszek  Demkowicz 
University  of  Texas  at  Austin 


The  report  describes  a  computational  and  theoretical  work  toward  developing  a  fully 
automatic  hp-adaptive  methodology  for  solving  3D  time-harmonic  Maxwell  equations. 
The  algorithm  produces  automatically  a  sequence  of  hp  meshes  that  deliver  exponential 
convergence.  The  optimal  hp-refinements  are  identified  by  minimizing  the  projection 
based  interpolation  error  of  the  solution  obtained  on  a  (fine)  grid  obtained  by  a  global  hp- 
refinement  of  a  current  (coarse)  grid. 

The  idea  of  the  algorithm  is  very  general  and  it  applies  to  all  discretizations  forming  the 
hp  de  Rham  diagram  -  HI-  H(curl)-and  H(div)-conforming. 


Final  report  for  the  project  on 


Fully  Automatic,  hp~ Adaptive 
Simulations  for  Maxwell’s  Equations. 

1  Accomplishments  of  the  project.  A  summary 

The  project  was  a  continuation  of  an  earlier,  three-year  project  on  ”  An  hp  Adaptive  Finite  Element 
Simulator  for  Maxwell’s  Equations  with  Application  to  3-D  Scattering  and  Resonating  Cavities 
Problems”,  during  which  a  new  Finite  Element  Method  allowing  for  variable  mesh  size  h  and 
order  of  approximation  p,  was  proposed  and  studied.  The  new  element  was  shown  to  belong  to  a 
general  family  of  conforming  elements  forming  an  exact  sequence,  and  a  family  of  projection-based 
interpolation  operators  was  constructed  that  make  the  de  Rham  diagram  commute.  Two  pilot 
codes,  2D  (triangles, quads),  and  3D  (hexahedra)  were  built,  and  a  number  of  model,  test  cases 
was  solved.  The  possibility  of  achieving  an  exponential  convergence  for  problems  with  singular 
solutions  was  demonstrated  through  numerical  experiments. 

In  the  presented  project,  the  crucial  task  of  working  out  a  fully  automatic  /ip-adaptive  strat¬ 
egy  that  would  deliver  the  exponential  convergence  for  solving  time-harmonic  Maxwell  equations, 
was  undertaken.  The  work  combined,  in  a  very  “tied  loop” ,  development  of  a  new,  fundamental 
mathematical  theory,  with  a  professional  quality  code  development.  The  two  sides  of  the  project  - 
theory  and  code  building,  strongly  influence  each  other.  Without  any  exaggeration,  it  would  have 
been  impossible  to  develop  the  projection-based  interpolation  theory  without  the  experience  gained 
dmring  the  computations,  with  the  theory  in  turn,  providing  theoretical  foundations  for  the  mesh 
optimization  algorithm  . 

1.1  Theory 

The  concept  of  elements  of  variable  order  forming  the  exact  sequence,  with  the  projection-based 
interpolation  operators  satisfying  the  commuting  diagram  property,  was  generalized  to  Nedelec 
simplices  of  the  first  kind  [4].  Out  of  all  Nedelec  elements,  only  one  (hexahedra  and  quads  of  the 
second  type),  do  not  satisfy  the  exact  sequence  property,  and  fail  to  converge.  Generalizations 
of  all  remaining  elements  to  the  variable  order  of  approximation  setting,  is  discussed  in  [11,  12]. 
The  optimal  p  and  hp  interpolation  estimates  were  derived  in  [5]  for  triangular  meshes.  The  3D 
interpolation  theory  is  much  more  advanced  and  different  principles  have  to  be  used  [13,  14].  While 
the  interpolation  theory  for  elliptic  problems  is  complete,  the  theory  for  H{curl)-  and  H{div)- 
conforming  discretizations  hinges  on  a  single  conjecture  on  existence  of  3D  polynomial  preserving 


extension  operators  for  H{curl)  spaces. 

Extensions  to  non-simplicial  elements  are  rather  straightforward,  and  we  plan  to  do  it  in  a  near 
future. 

In  a  separate  effort,  we  have  been  studying  the  discrete  compactness  property  for  the  p  and 
hp  methods.  This  crucial  property  implies  convergence  of  eigenvalues  and,  therefore,  stability 
of  FE  discretizations  of  time-harmonic  Maxwell  equations  as  well.  Some  preliminary  results  for 
a  triangular  element,  obtained  under  a  -discrete  stability  conjecture,  were  reported  in  [2].  I 
am  happy  to  inform  that  the  results  have  been  recently  extended  to  the  quad  element  (paper  in 
completion) . 

1.2  Code  development  and  problems  solved 

The  original,  proof  of  concept  codes,  developed  with  Rachowicz  [3],  have  been  rewritten  within 
a  completely  new  data  structure  [7,  8].  The  codes  survived  the  ultimate  test  of  supporting  the 
automatic  /ip-adaptivity  and  a  two-grid  solver.  Since  last  Summer,  the  codes  are  maintained  on 
the  Web  under  CVS  with  current  updates  available  immediately  to  computational  math  community. 
The  number  of  national  and  international  users  is  steadily  growing.  Most  of  the  coding  is  done  in 
collaboration  with  Waldek  Rachowicz  from  Cracow  and  Adam  Zdunek  from  Stockholm. 

The  two  principal  components  of  this  project  are  the  automatic  /ip-adaptivity  module  and  the 
two  grid  solver.  The  2D  mesh  optimization  module  has  undergone  four  iterations  [6,  9,  10]  and  it  is 
operational  for  both  elliptic  and  Maxwell  problems.  The  methodology  has  been  successfully  tested 
on  tvelwe  non-trivial  test  examples  including  orthotropic  heat  conduction,  elasticity,  acoustic  and 
electromagnetic  scattering  problems,  and  wave  propagation  in  conductive  media. 

It  took  two  and  a  half  years  to  complete  the  3D  mesh  optimization  module  which  has  been  so 
far  tested  only  on  two  elliptic  examples.  A  second  implementation  of  the  3D  mesh  optimization 
module  is  underway  (a  Ph.D.  project  of  Jason  Kurtz).  We  anticipate  presenting  solutions  to  the 
3D  scattering  problems  obtained  with  the  automatic  hp-refinements  algorithm,  during  the  next  Air 
Force  workshop  (2005). 

Crucial  to  the  success  of  the  proposed  methodology  is  the  two-grid  solver  that  makes  the  solution 
of  the  fine  grid  problem  possible.  The  work  has  been  summarized  in  [17]  and  in  dissertation  of 
David  Pardo  [18],  to  be  defended  on  April  2,  2004. 

In  summary,  on  the  code  development  side,  we  are  about  a  year  late,  compared  with  the  initial 
plans. 

In  a  parallel  effort,  the  automatic  /ip- adaptivity  algorithm  has  been  coupled  with  the  goal- 
oriented  mesh  optimization  [20],  and  the  2D  elliptic  code  has  been  successfully  parallelized  for 
memory  distributed  platforms. 
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